Pattern Formation in a Simple Model of Photochemical Reaction 
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Stripe and columnar patterns were experimentally found in photopolymer films during irradia- 
tion of collimated UV light. We propose a mathematical model of photopolymerization, perform 
numerical simulation, and find a spatially periodic columnar structure. 
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Recently, the self-organization of microstructures has been intensively studied in physics, material science, and 
chemistry P, Q Self-organization and pattern formation have also been studied in polymer science. There are various 
mechanisms of self-organization in polymer science. Spatially periodic Turing patterns were found in reaction-diffusion 
systems in gels. 0, |j] Honeycomb structures were formed when an amphiphilic polymer solution was cast on solid 
surfaces, which was due to the Marangoni convection. [5j Various patterns were created on the surface of polymer 
gels during the volume phase transition owing to mechanical instability. [6| Various mesoscopic patterns were found in 
block copolymers by micro phase separation. Q 

In this study, we consider self-organization in photopolymerization. In photopolymerization, UV light is illuminated 
onto a monomer solvent including photoinitiators. The initiators are excited by UV light and radicals are created. A 
monomer with a radical is added to a polymer and the polymerization chain grows. Photopolymerization stops when 
radicals become inactive or markedly decrease in number. There are many applications of photopolymers, i.e., a UV 
curing resin and a photoresist for lithography. Various types of pattern formation were reported in photo-crosslinking 
and photopolymerization reactions. For example, Nakanishi et al. found a hexagonal phase during the synthesis 
of interpenetrating polymer networks in a polymer mixture. They concluded that pattern formation occurs by the 
cooperation between phase separation and elastic repulsion. [8| 

On the other hand, Honda et al. found stripe patterns in photopolymer films during the irradiation of tube UV 
light, m Uematsu et al. found amorphous columnar patterns using point light. flol | In their experiment, a monomer 
mixture of methacrylate monomer and photoinitiator was injected into the space between two coverglasses, and UV 
light was irradiated onto the setup. The thickness of the space was controlled. No patterns appeared in a thin 
sample of 0.45 mm, but patterns appeared in a thick sample of 1.5 mm thickness. In the thick sample, a pattern 
was generated below a depth of about 200 /xm from the irradiated surface. The wavelength of the stripe pattern is 
about 10 [im. Chatani and Hattori performed a similar type of experiment using different monomer mixtures and 
confirmed columnar patterns. They found that no pattern appeared in a thin sample and a clear pattern in a medium- 
thickness sample; however the pattern became random in a thick sample. [IlT | Pattern formation seems to occur during 
photopolymerization but its detailed mechanism is not yet understood. Honda et al. suggested a diffraction effect 
and Uematsu et al. proposed a type of lens effect as a mechanism of pattern formation. That is, if the intensity 
of UV light is inhomogeneous as a result of some perturbations, photopolymerization proceeds inhomogeneously and 
the spatial inhomogeneity of the refraction index increases, which makes the inhomogeneity of the UV light intensity 
increase further by the diffraction and (or) refraction effect. This type of instability in the photochemical reaction 
has not yet been studied mathematically. 

We numerically study a mathematical model of photochemical polymerization. We assume radical polymerization 
induced by photoinitiators. Firstly, radicals, R, are generated by the decomposition of initiators, I, by UV radiation. 
The addition of a radical to a monomer, M, produces a propagating radical, Q. The chain propagation occurs by the 
addition of a monomer to the propagating radical. The reaction of two propagating radicals terminates the chain 
propagation and yields a polymer, P. The reactions are expressed as 



We assume that the rate of the production of radicals from initiators is proportional to the light intensity \4>\ 2 . The 
reaction rate equations are expressed as 



I -)• 2R, 



R + M -> Q, 
Q + M Q, 
2Q -> P 



(1) 




2k 4 Q 2 , 



k 2 RM, 



k 3 MQ, 



(2) 



where I , R, ■ ■ ■ , P are the concentratios of I, R, - ■ •, and P, and k\, ■ ■ ■ , k^ are reaction rate constants. The concentration 
of I decreases monotonically and becomes zero finally. The polymer concentration P increases monotonically and 
becomes constant by the termination of the reactions. The stationary state depends on the initial conditions, which 
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FIG. 1. (a) Polymer concentration P as function of |<j!>| 2 in eq. (1) for the three initial conditions of 1(0) = 0.1,0.3, and 
1. The other initial conditions are M(0) = 1 and i?(0) = Q(0) = -P(O) = 0. The parameter values are ki = 0.005, ki = 
kz = ki = 1, and Di — Dr — Dm = 1. (b) Amplitude W of spatially periodic pattern obtained from the initial condition 
I(x,z) = 0.998 — 0.002 cos(fca;) with various wavenumbers k. (c) Pattern of polymerization. In the shaded region, P > 0.2625. 
The initial condition of I is random between 0.9 < I < 1. 



is rather different from the simple model studied in previous sections. This is because the polymerization reaction 
expressed by eq. (1) is irreversible. 

We have performed numerical simulation of eq. (1). Figure 1(a) shows P in the stationary state for k\ = 0.005 
and k 2 — k 3 — fc 4 = 1 as a function of the light intensity \4>\ 2 ■ The initial conditions of I are set to be 7(0) = 1, 0.3, 
and 0.1. The other initial conditions are M(0) = 1, R{0) = 0, Q(0) = 0, and P(0) = 0. The stationary value of P 
increases rapidly near \(j>\ 2 = with the light intensity for the initial conditions, because the production of radicals 
from photoinitiators is proportional to the light intensity. 

To study pattern formation, we investigate the following coupled reaction-diffusion equations: 



m 

dt 
dR 

~dt 
dM 

~dt 
dQ 
dt 
dP 
~dt 



-ki\4\ 2 I + DiV 2 I, 
2k 1 \(j)\ 2 I - k 2 RM + D R V 2 R, 
-k 2 RM - k 3 MQ + D M V 2 Q, 
k 2 RM ~ 2fc 4 Q 2 , 

fc 4 Q 2 , (3) 



where Dj, D R , and D M are diffusion constants for I, R, and M, respectively, and the diffusion of macromolecules Q 
and P is neglected. These reaction-diffusion equations are coupled with the Schrodinger-type equation for the light 
propagation: 

4^-^-aP0. (4) 
oz 2 

Here, we have assumed that the refractive index is proportional to the polymer concentration P. 

We have performed numerical simulation in a two dimensional system of size L x x L z . The parameters are set 
to be k\ = 0.005, k 2 = k$ = fc 4 = 1, Dj = Dr = Dm = 1, and a = 200. The input light intensity is assumed 
to be uniform (f>(x, z) = 1 at z = 0. A stationary pattern of P is obtained from the initial condition of I(x, z) = 
l-0.1{l+cos(247rx/L K )}, M(x, z) = 1, and R(x, z) — Q(x, z) — P(x, z) = for L x xL z = 80 x 15. A spatially periodic 
pattern with k = 24w/L x is created from these initial conditions. We have calculated the amplitude W of spatially 
periodic structures between 0<z<z = 5byW= [f*° J Q m P 2 dxdz/(z L x ) — {f*° J Q x Pdxdz) / (zqLx)} 2 ] 1 / 2 for 
the initial conditions I(x,z) = 1 — 0.002{1 + cos(fcx)}, M (x, z) = l,R(x,z) = Q(x,z) — P(x,z) — with various 
initial wavenumbers k. Figure 1(b) shows the dependence of W on k. The amplitude W takes its maximum near 
k = 22ir/L x . This implies that a pattern with a finite wavenumber near k = 22tt/L x grows dominantly. 

We have performed numerical simulation from a random initial condition. That is, I(x, z) is assumed to be a random 
number between 0.9 and 1, and the other initial condition is the same as before: M(x,z) = l,R(x,z) = Q(x,z) — 
P(x,z) — 0. The input light is uniform in this numerical simulation, that is, (/)q(x) = 1 at z = 0. Figure 1(c) 
shows a stationary pattern of P for L x x L z = 80 x 15. A spatially periodic columnar structure of about 12 columns 
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FIG. 2. (a) Pattern of light intensity in thick system of L x x L z — 40 x 60. In the shaded region, \4>\ 2 > 0.2. The initial 
condition of / is 1, and the input light is (f>o(x) = 0.3/cosh{(x — L x /2)/6}. The parameter value of a is 80. (b) Pattern of 
polymer. P > 0.6 in the shaded region. 



appears, although some randomness is overlapped. The wavenumber is therefore k ~ 24-7T j L x . The spatially periodic 
columnar structure is clearly seen only for z > 2. Even if the system size is twice enlarged to L x x L z = 160 x 15, 
we have obtained a similar pattern of 24 columns. We think that the fact that a pattern with a finite wavenumber 
with k r~j 2Air/L x is self-organized from a random initial condition is closely related to the numerical result shown in 
Fig. 1(b) that a pattern with a finite wavenumber tends to grow fast. 

We have performed another type of simulation for the self-organization of a light waveguide in this thicker system. 
In this simulation, the input light is not uniform, but localized in space. We have assumed that the input light is 
localized as 4>o(x) = 0.3/cosh{(a; — L x /2)/6}. Figures 2(a) and 2(b) respectively show stationary patterns of |</>| 2 
and P, where \<fi\ 2 > 0.2 and P > 0.6 are satisfied in the shaded regions. The parameter a is 80 in this numerical 
simulation. The beam structure and strongly polymerized region are localized near x = L x /2 and rather uniform in 
the z-direction. This localized pattern corresponds to a soliton state. 

Similar type of self-organized light waveguide was experimentally studied in the past by Kagami et al., who found 
that a light waveguide self-organizes by photopolymerization with beam illumination. [12j We think that their self- 
organized light waveguide is closely related to our simulation result. 

In summary, we have proposed a model of photopolymerization. The irreversibility of photopolymerization is 
incorporated in this model. The reaction rate of the decomposition of initiators is assumed to be proportional to 
the light intensity, and the complicated reaction-diffusion equations are coupled with the Schrodinger-type equation 
for light propagation. The complicated model equations are difficult to analyze theoretically, but we have performed 
numerical simulation. A spatially periodic columnar pattern is obtained from a slightly random initial condition of 
I(x, z), which is qualitatively similar to the pattern formation in the simple model. 
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